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Abstract 

Some problems related to the transition density u(t, x) of the diffusion on the 



^ ■ Sierpinski gasket are considerd, based on recent rigorous results and detailed 

^^ ■ numerical calculations. The main contents are an extension of Flory's formula 

for the end-to-end distance exponent of self-avoiding walks on the fractal 
G^ ■ spaces, and an evidence of the oscillatory behavior of u{t, x) on the Sierpinski 



gasket. 
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In this letter we report our study on the transition density u{t,x) of the diffusion and 
the random walk on the Sierpinski gasket, based on recent rigorous results and detailed 
numerical calculations. The main contents are an extension of Flory's formula for the end- 
to-end distance exponent of self-avoiding walks on the fractal spaces, and an evidence of the 
oscillatory behavior of u{t,x) on the Sierpinski gasket. 

Recently, rigorous justification of the (symmetric and isotropic) diffusion on the Sierpin- 
ski gasket and analysis of its behavior appeared in mathematics literatures |l]J^. Among 
other results in these studies, we focus on the transition density u(t, x), the density at point 
X at time t > 0, for the diffusion starting at t = from the origin of the Sierpinski gasket. 
In [0 u{t, x) is rigorously shown to exist, and the following form of bound is proved to hold 
for all t > and at any point x on the Sierpinski gasket: 

fit,x;C,,C2) < uit,x) < fit,x;C3,C,) , (1) 

where CjS are some positive constants independent of t and x, and the function / is given 
by 

fit, X- Ci, C2) = Cir^^/2 exp{-C2(|a;|ri/'^-)2''} , r^ = d^/{2d^ - 2) . (2) 

The exponents d^j and dg in (^ are the walk dimension and the spectral dimension, re- 
spectively, whose values for the Sierpinski gasket are d^j = log 5/ log 2, and dg = 2df/dw = 
2 log 3/ log 5, where dj = log 3/ log 2 is the fractal dimension @,^. The specific form \x\t~^^'^^ 
in (^ implies anomalous diffusion < |a;(t)p >~ t'^/'^^" (or rather, this relation defines d^). 
Note that the value of 77 in (Q) cannot be determined from this relation alone. Bounds of the 
form (|l|) with 77 as in (^) are mathematically proved to hold also for a wide class of finitely 
ramified fractals ||5[ (with some generalizations which we will not deal with here), and even 
on some infinitely ramified fractals such as the Sierpinski carpet |^. The wide applicability 
of (0) suggests us to take this formula as one of the basis in the studies of u(t,x). We 
consider two problems related to u{t,x). One is the extension of Flory's formula for the 
self-avoiding walks (SAW) to the fractal spaces, and the other is the oscillatory behavior of 
u(t,x) on the Sierpinski gasket. 



Consider a SAW on a fractal with the fractal dimension df and the spectral dimension 
dg. The end-to-end distance exponet u is defined by R{N) ~ N'^ {N ^1), where A^ is the 
number of steps of a SAW and -R(A^) =< |a;(A^)| > is the average end-to-end distance of 
A^-step SAW. According to the mean-field type arguments for SAW, Flory's value up 0)111 fo^ 
the exponet u is obtained by finding the solution R = Rp{N) which attains the minimum 
of the 'free energy' —\ogu{N,R) + V{N,R) for each A^, where we wrote u{N,R) for an 
average of the transition density u{N,x) of the simple random walk over x with |x| ~ R, 
and V^(A^, R) = N"^ / R'^f represents the volume exclusion effects, vp is then determined by 
Rf{N) ~ A^*^^. The studies that derived (0) for finitely ramified fractals start with analysis 
of simple random walks and then reach the diffusions by taking continuum limits. Therefore 
the long time behavior (A^ ^ 1) of the transition density u{N,x) for a random walk also 
satisfies (P with (Q). We use the form (0) for u{N,R) to obtain, 

^F = M^) = 2 ^/^{f" , r] = dJ{2d^-2). (3) 

df + 2r] 

The argument holds for any network with definite fractal dimension df and walk dimension 
dw The value r] = df/dg = dy^/2 was proposed at times when the form of (^) was not 
settled, resulting in a simpler formula h'F{dw/2) = 3/{df + d^) HJ^- In [^ it was pointed 
out that there was no justification in this choice other than simplicity, and that the problem 
of the choice of r] remained open. A heuristic explanation of the rigorous proof p| for the 
value f] = dw / {2dyj — 2) is as follows. 

Fix the step A^ and the distance R such that R= \x\ ^ N^^^^' and consider walks of A^ 
steps that reach a point at distance R; i.e. look at walks going outwards quickly. Classify 
the random walk sample paths by the scale tq, such that for scale r larger (resp., smaller) 
than To the random walker walks straight (resp., walks randomly; AA^ ~ Ar'^™). A walk 
specified by the scale tq passes straight through R/r^ blocks of scale tq, by definition of tq. 
Since it takes steps of order rg™ to pass through each block, we have A^ ~ -R/tq ■ Tq"" , which 
imply that the dominant contribution to the quick diffusion specified by (A^, R) comes from 
the walks with tq ~ (N/RY^^'^^~^\ Each time the walker passes the block straight through 



he loses probability by 1/4, because, at each node, there are 4 possible directions (i.e. the 
4 outmost vertices of the two blocks connected to the node) to go. The total decay of 
probability, which gives an estimate of the transition density is u{N, R) ~ 4"^^/'"''^ because 
the walker passes R/ro blocks straight through. Using the estimate for r^ given above, we 
have — logpAr(-R) ~ (^Ju^-'^-fdwyw/idn,-!) log 4, which implies t] = d^/{2dw — 2). 

The reason that (0) is to be used for Flory's formula can be seen from the above argument: 
SAW is 'pushed outwards' compared to random walks, owing to self-repulsion or volume 
exclusion effects. Therefore, the dominant contribution to SAW comes from those walks that 
move quickly away. The argument given above explains that the value rj = dw/{2duj — 2) 
is the consequence of the contribution from walks which quickly move away, hence it is 
reasonable to use this form in deriving the Flory's formula. 

The explicit values of df and d^ are known for the Sierpinski gasket and its natural 

(i-dimensional generalizations (dSG), constructed by {d + l)-simplex instead of triangle to 

construct the Sierpinski gasket. The values are dj = log{d + l)/log2 and d^ = \og{d + 

3)/ log 2 for dSG ||5,|], with which up = upi^dyj / {2dyj—2)) can be calculated from (|^). For 2SG 

(= Sierpinski gasket) and 3SG, the values of up (0.8249 • ■ ■ and 0.724588 ■ ■ ■, respectively) 

are to be compared with the exact values of z/, which are 

v{2SG) = log 2/ log (7 - v^)/2 = 0.79862 ■ • • , 

(4) 
T^{3SG) = 0.67402 • ■ ■ . 

The value of ulSSG) has similar exact expression as that for i'{2SG), but in place of integers 

7 and 5 appear roots of a 14-th order algebraic equation |1^ . The values in (H) have been 



known for some time |jTT[ (see also |]8|JT^), and are recently proved rigorously in [|T3|,|T0 



Flory's formula (^ is within 3% and 8% precision from the exact values for 2SG and 3SG, 
respectively. Flory's formula is known to be numerically very good for SAW on Euclidean 



lattices (for a recent nice review of SAW on Euclidean lattices, see IjT^)- The extended 
Flory's formula (^ which we have is not very bad, but not very close to the exact values 
compared to Euclidean cases. If we put r] = 1, the values become closer (in fact, it is close 
to the best choice) to the exact results (0.06% and 3% deviation for 2SG and 3SG, resp.). 



The choice rj = dw/{2dw — 2) has the most sound basis (HD and (0), but the value is better 
for 2SG and 3SG with rj = 1. Deviation of (^ from the rigorous and exact results (|) leads 
our interest to detailed numerical studies oiu{t,x). 

Two open problems are found in the literatures concerning the detailed structure of 
u{t,x). One problem is the value of 77; the results given in |]IB| agree with (^, while those 



in |]T6| claim the value rj = dy^/2. The other problem is the observation in [@,0 that there 
are 'oscillations' in u{t,x). 

To perform numerical calculations, we regard u{t,x) as the electric charge density of a 
point X at time t, and reformulate the problem in terms of the impedance circuits. We 
consider a d- dimensional Sierpinski gakset dSG. The corresponding electrical circuit has 
impedance distribution on the gasket, and also impedance between the gasket and the 
ground. By symmetry and star-triangle {Y — A) type relations, a unit block of the gasket (a 
sub-circuit of (i-dimensional simplex with side length 1) can effectively be represented by a 
device with rf + 1 terminals, each connected by an impedance a(s) to the center point of the 
simplex, to which the ground is connected by an impedance &(s). s is the dual variable to 
t in the Laplace transform. The self-similarity of the diffusion implies the scaling behavior 
B u{Lt, 2x) = u{t, x) [^^, where we put L = d + 3 = 2'^^ and B = d + 1 = 2'^i . Using the 
scaling behavior, the self-similarity of the gasket, and the similarity among d + 1 terminals 
of a block simplex, together with star-triangle type relations, we find 

{a{Ls) L/B, g{Ls)) = W{a{s), g{s)) , (5) 

where g{s) = 2a{s)/b{s), and W{x,y) = {x{y + L)/{y + B),y{y + L)). We also find g{0) = 0. 
g'{0) and a(0) determine the normalization of u and t. We focus on the normalization 
independent quantities such as exponents and oscillations. By fixing g'{0) and a(0), the 
solution of the functional equation (^ is determined uniquely. The equation for g{s) is 
known as Schroder's functional equation. The existence of the solution has been studied 
|T^, but its detailed behavior seems to be unknown. 

We define two asymptotic functions C{s) = g^^'^s/^ \[tq^ (L/B)"'a{L"'s) and k{s) = 



s ^''^"' lim 2 " log g{L'^s), for s > 0. These functions are periodic in logs with period 
logL, hence can be expanded in Fourier series: 

oo 

C{s) = Co + XI Cn sin(27rn log^ s + (j)n) , 

'Si (6) 

k{s) = ko + ^kn sin(27™ log^ s + (f)'J . 

n=l 

We numerically obtained by double precision FORTRAN calculations for d = 2, the 
Sierpisnki gasket; ci/cq = 1.21929438 x 10-^ and ca/ci = 3.68 x 10"^ for C{s), and 
ki/ko = 1.5264191 x 10~®, and ^2/^1 = 5.6 x 10~^, for k{s). The results show a strong hi- 
erarchy of coefficients, such that the higher frequency components have exponentially small 
values {cn = O{10~^^), kn = O(10~^")). Note that we have small but non-zero numbers of 
0(10^^) out of equation (|) with 0(1) coefficients, which is potentially an interesting phe- 
nomena. We performed the numerical calculations up to d = 10, and obtained qualitatively 
similar behavior, with somewhat larger amplitudes of oscillations for larger d. For precision 
check, we performed the calculations for d = 1, corresponding to the diffusion on a line, and 
obtained the correct constant values C{s) = 2^/^ and k{s) = 2^/^, within error 10"^^. The 
peak values of k{s) are consistent with the numerically obtained values in 0. 

We can calculate the Laplace transform u{s, x) of the density using the impedances. For 
the transition density at the origin, we have u{s,0) = aoo(s)/2, where aoo(s) = C{s) s'^'^/^-i 
is the impedance corresponding to a simplex of 'infinite' size. Using (P), we can evaluate the 
inverse Laplace transformation of u{s, 0), term by term in the Fourier series, using change 
of contours. We have 

00 

U{t,0) = Nt~'^''/'^[CoPo + Y^ CnPnCOs{27m\ogLt - iJn)] , 

n=l 

with Cn as in (D and p„ = \T{2-^ds + 7r-^nni)\2-^/^{cosh{2nn) - cos(7rrf^))i/2 for n > 0, 
where Q = 27r^/logL. A^ is a normalization constant. Thus the oscillation in C{s) explains 
that in u{t, 0). 

We parametrize u{t,x) = f(t,x;Ci(t),C2(t,x)) with / as in (|^), and consider the 
oscillations in CjS. Using the values given below @ for c„s, we have, for d = 2, 
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Ci{t) = t'^^/'^u{t, 0) = Cio + Cii cos(27r log^ t - ^i) + ■ ■ -, with Cu/Cw = 8.0964779 x IQ-^. 
We performed numerical Laplace inverse transformation of u{s, 0) and obtained a consistent 
value. The Laplace transform of the density u{t, 1) at a vertex of a unit simplex is given 
by u{s,l)/u{s,0) = (1 + gis)/i2d)){l - a{s)/ao.{s)) - a{s)/{da^{s)). C2{t) = C2(t, 1) is 
then given by C2{t) = t^''/'^™ log(M(t, 0)/u{t, 1)). To calculate the inverse Laplace transforms 
for small t {t < 10^^), we find s = Sq > which gives minimum of u{s, l)exp(st) (as we 
do in the steepest descent method of complex contour integration), and numerically eval- 
uate the contour integration with the contour 9fJ(s) = Sq. For larger t we use the contour 
3fJ(s) = 1/t. (The details of the numerical calculations will be reported elsewhere.) Figs. [1| 
and |] show C2{t) for (i = 2 as a function of Int. Fig. |^ shows the small oscillation (O(10~^) 
in amplitude) in C2{t) for very small t. We performed similar calculations for the exactly 
solvable d = 1 case, and checked that the error in the range of Fig. |^ is 0(10^^^). Small 
t corresponds to large s in the Laplace transform, where we have an asymptotic formula 
u(s, l)/u(s, 0) ~ exp(— S"'^/'^™/i:(s)). Thus the oscillation in k{s) explains that in C2{t) for 
very small t. 

From Fig. |l] we see that (besides the tiny oscillation) C2{t) is flat for small t. This is 
consistent with the fact that the value rj = du,/{2dw — 2) explains the asymptotic behavior 
of u{t,x) as t ^ 0. For larger t, C2{t) is decreasing with a significant size of oscillation. If 
we try to explain this decrease in terms of 'effective (dynamical)' changes in the value of 1], 
i.e. keep C2 constant and let t] change as 77 = rjeffit), we see that the effective value rjeffif) 
increases as t is increased. Our data are in favor of the argument in |TB[ that f]eff{t) = d^/2 



T6| may be effectively good for t = 0(1), while r] = dw/{2dw — 2) is good for t <^ 1. Note 
also that du,/{2d^ — 2) < 1 < dw/2, where rj = 1 gives good i^f{v)- The oscillatory behavior 
of the data prevents us from obtaining precise results on the value of ^e//- In contrast to the 
clarity in the meaning of the value rj = dw/{2dw — 2), theoretical basis for r/g// is still unclear, 
which has to be settled before we can be conclusive about its value and implications. 
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FIGURES 
FIG. 1. C2{t) = t^Vrf^ log{u{t,0)/u{t, 1)). 

FIG. 2. Oscillation of C2{t) for small t. 
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Fig. 1 
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